SPT-03/139 



The PDEs of biorthogonal polynomials 
arising in the two-matrix model ^ 



O . M. Bertola^^^ B. Eynard^'^^ 

^ I ' Centre de recherches mathematiques, Universite de Montreal 

' CP. 6128, succ. centre ville, Montreal, Quebec, Canada H3C 3J7 

^ Department of Mathematics and Statistics, Concordia University 

7141 Sherbrooke W., Montreal, Quebec, Canada H^B 1R6 

* Service de Physique Theorique, CEA/Saclay 
^ ■ Orme des Merisiers F-91191 Gif-sur-Yvette Cedex, FRANCE 



(N 

> . 

, Abstract: 
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, The two- matrix model can be solved by introducing bi-orthogonal polynomials. In the case the potentials in the 
measure are polynomials, finite sequences of bi-orthogonal polynomials (called windows) satisfy polynomial ODEs as 
' well as deformation equations (PDEs) and finite difference equations (AE) which are all Frobenius compatible and 
" define discrete and continuous isomonodromic deformations for the irregular ODE, as shown in previous works of ours. 
In the one matrix model an explicit and concise expression for the coefficients of these systems is known and it 
allows to relate the partition function with the isomonodromic tau-function of the overdetermined system. Here, we 
^ [ provide the generalization of those expressions to the case of bi-orthogonal polynomials, which enables us to compute 
' the determinant of the fundamental solution of the overdetermined system of ODE+PDEs+AE. 
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1 Introduction 



The two-matrix model and biorthogonal polynomials have recently witnessed a renewed interest due to the hope that 
the asymptotics could be found from a Riemann-Hilbert approach similar to that used in one-matrix models ([^1 
|12[I13| . |17| . |S] 0], Uni)- Our recent papers have shown some unexpectedly rich algebraic structure associated to the 
biorthogonal polynomials for two polynomial potentials [51 0] . We recall that the two matrix- model is defined as 
the ensemble of pairs of hermitian matrices of size N, with the measure: 

^g-i[tr Vi(Afi) + V.(Af.)-MiM2]dMidM2 (1.1) 

where dMidAf2 is the standard Lebesgue measure (product of Lebesgue measures of all real components of Mi and 
M2), V\ and V2 are two polynomials (called the potentials) of degree di + 1 and di + 1: 

^i(-)-Et^'' ^^(^)=Ef^'= (1-2) 

fe=l fe=l 

and 2,]^ is the partition function: 



The two matrix model has a wide range of applications in physics, in particular Euclidean 2d gravity i.e. the statistical 
physics of a random surface |^ El ■ 

The 2-matrix model can be solved with the help of two families of bi-orthogonal polynomials 

7r„(a;) + , cr„(y) = + . . . (1.4) 

which are such that: 

y y dxdy ^„(x)a™(y) e-*[^i(-)+^^(«)--^l = /i„(S„™ . (1.5) 

It was shown in that a sequence of (i2 + 1 consecutive tt^'s (or di -I- 1 consecutive cr^'s) obeys a closed ODE 
with respect to the "spectral parameter" x (or y) 



h— 

da; 



Di{x) 

N 



\ '>Pn{x) 



(1.6) 



as well as deformation equation (PDEs) with respect to infinitesimal variations of the coefficients of the potentials 
Vi , V2 and difference equations with respect to the choice of the consecutive elements ( AE) . The aim of this article is 
to provide for an explicit expression of said ODE and PDEs, which is similar in spirit to the work done in ll()| (the 
AE are already very simple and can be found in |5j ) . 

Here we follow the same logical lines that were followed in our previous 0. 

We should remark that an explicit expressions for the ODE is given in 5 but is not convenient for practical 
purposes. For instance it is not at all obvious how to compute the trace of the matrix Di{x) from the expression 

N 

written in |5l . We expect that all the spectral invariants of Di (x) will play an essential role in the future developments 
on the basis that in the one-matrix model (HllHl the corresponding invariants allow to put in direct relation the partition 
function of the model with the isomonodromic tau- function in the sense of ^1 ^1 El • Indeed it is well known that 
the one-matrix model is a subcase of the two-matrix model where one of the potential is Gaussian. 

The expressions contained in this paper allow us to easily compute the first invariants (traces) of the ODE's and 
PDE's, thus allowing us to evaluate the determinant of a fundamental system of solutions of the Frobenius compatible 
system. Hopefully they will prove useful for the higher invariants as well. 
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Remark 1.1 (Notation) In the following we will often consider semi-infinite matrices. If M is a semi-infinite 
matrix we denote: 

• Af+ the strictly upper triangular part of M , 

• M_ the strictly lower triangular part of M , 

• Mk the fc*'' diagonal above the main diagonal if k>0, and the diagonal below the main diagonal if k <0, 

• Af+o = M+ + iMo, 

• M_o = M- + ^Mo, 

• Af> the upper triangular part (including the main diagonal), 

• Af< the lower triangular part (including the main diagonal). 

We introduce the following semi-infinite matrices: 
A is the shift matrix: 

Ay = Si+ij . (1.7) 
n is the projector on the subspace spanned by the first + 1 basis vectors: 

N 

N + l-times 



n = diag( 1,1,..., 1 ,0,...) . (1.8) 



N 

We have the simple properties (L has only a right inverse): 

An= n A , nA* = A* n (1.9) 

N N-l N N-1 

A*A = l-n, AA* = 1 , An = o, nA* = o. (i.io) 

00 

Finally the notation for the folded matrices is slightly changed from ^1 in order to make it more consistent ( in 
particular the location of certain over /under- scripts). 

2 Biorthogonal polynomials: definitions and properties 

It is convenient to introduce the wave functions and their Fourier-Laplace transforms, related to the bi-orthogonal 
polynomials H1.5() by the relations: 

^^{x) - ^7r„(x)e-*^i(-) , My) = ^a„(y)e-*^=(«) (2.11) 
±Jy) = J dxe^-yM^) , ^J^) = J d2/e*->„(y) . (2.12) 

The unspecified contour of integration can be chosen in di homologically distinct ways for the functions V'„(y) ('^2 
ways for (j>^{x)) as explained in j^Ej. We introduce also the semi-infinite wave- vectors 

vI/ = (V.o>i,...)' , a> = (00,01,...)* , * = (V'o,V'i,...) , $ = (00,01,...). (2.13) 

oo oo oo oo ~ 

In order to express the ODEs that these wave- functions satisfy we define the following "windows" of consecutive 
wave-functions of dimension ^2 + 1 (resp. di + 1): 

5' = (V'AT-da, . . . , V'Ar-i,'0Af)* , $ = (0JV-di,...,0Af-l,0Ar)* (2.14) 
N N 



N N 



-N-V2-N 



7 V Ar J • ■ ■ I ' 



(2.15) 



The two windows in each of the pairs ^1/,$ and $, j!. (each constituted of windows of the same size) are called 

\ff J \N J 

"dual windows" as explained in |S] . It is shown in |S1 0] that -similarly to the ordinary orthogonal polynomials- 
the semi-infinite wave- vectors obey: 



oo oo oo oo oo oo 



(2.16) 



oo oo 



?i9^$ = p^, - hdy <P Q* <i> , hdy^^^Q , hdy^ = ^P* (2.17) 

CxD oo oo OO OO 



OO oo 



where the semi-infinite matrices P and Q are of finite band size. Componentwise the above equations translate into 

xipnix) = 7(n)V'«+i +'^akin)il)n-kix) , y(t>„{y) = j{n)(f>„+i f3k{n)(f>n^k{y) (2-18) 



A:=0 



fe=0 



d2 



x0^(a;)=7(n-l)0^_^+^afc(n + /c)0 (x) , #„(?/) = + /^'^('^ + %„+fc(y) (2-19) 



fe=0 



where ^{n) - 

Since they implement the multiplication and differentiation by the spectral parameters x or y, the matrices P and 
Q obey the Heisenberg relation: 



and they have the properties [H El ^ : 



[p, Q] = hi 

P+ = Vl{Q)+ , Q- = V^{P)- 
Po = Vl{Q)o , Qo = V^{P)o 

°n,„-i = i{n - 1) = F/(g)„,„_i - h 

Qn.n+l = 7(") = V2iP)n.n+l - tl 



l{n- 1) 
n + 1 



7(n) 



(2.20) 

(2.21) 
(2.22) 
(2.23) 

(2.24) 



The relations H2.17|l can be "folded" "F onto the finite windows defined in H2.14I2.15|I using the recursion relations 
(|2.16|l so as to give finite dimensional ODEs: 

N N N N N N 

-hd^^^ = Di{x)^^ , - hdy ^ = D2{y) ^ , fidy:^ = :^D2iy) , ;i9^$ = $^i(x) (2.25) 

N N N N N N 

N N 

where Di{x) and D.i{x) (resp. D2{y) and Daiv)) ^^'c matrices of dimension (i2 + 1 (resp. di -I- 1) whose entries are 

N N 

polynomials in x (resp. y) of degree at most di (resp. ^2)- They enjoy the duality property (l5 , Thm 4.2): 



N N N 



D^ix)A = ADi{x) , D^iy)] 

N 

N N 

where A and B are the Christoffel-Darboux pairing matrices: 



N N N 



( 



N 







ad,{N)... 

\ ad,{N + d2-l) 



N 



-7(iV-l)\ 
' 



/ 



(2.26) 



(2.27) 
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N 



( 







V l3dAN + di-\) 



'l{N-l)\ 
* 



(2.28) 



. N N 

The two C-D pairing matrices A and B arc the only non- vanishing blocks in the semi-infinite matrix commutators: 



[Q, n ] (resp. [P\ n ]) 



(2.29) 



where n is the projector defined in p. 8(1 . These matrices play an essential role in the computation of the relevant 



N-l 



spectral statistics of the model inasmuch as they enter directly in the kernels used in the Dyson-like formulas [T^ 



N N 

{x -x) Kii{x,x') = $(a;') A*(a;) = $(2:') 

N 00 
N N 

{y' - y) K22{y, y') = My') B ^(y) - *(y') 

N 00 



Q, n 
Q, n 

N-l 



(2.30) 



2.1 Deformation equations (PDE's) 



The deformations of the biorthogonal polynomials (derivations w.r.t the coefficients of the two potentials), are repre- 
sented by semi-infinite finite-band matrices: 



(2.31) 



00 00 



00 00 



0000 0000 00 cx; 0000 

where the semi-infinite matrices Uk, Vj are related to the matrices P, Q by the equations 



1 



which can be folded onto the windows using ea. H2.16|l as explained in |5] 



N N N N N N 

N N N ' N N N 



N 



(2.32) 



(2.33) 



(2.34) 



N N N N 1 _ _ 

hdu^ * = ^UK{y) , hd^j * = - * Vj(y) , hdu,, $ = ~UK{y) $ , hd,, $ = Vj{y) $ (2.35) 

N N N N N N 



where Uk{x) and U_x{'^) ^'^'^ matrices of dimension d2 + 1 whose entries are polynomials in x of degree at most k, 



N 



N 



and Vj{x) and V_j{x) are matrices of dimension c?2 + 1 whose entries are polynomials in x of degree at most Jdi (and 



N 



similar statement for the tilde-matrices with di replacing 6,2 and vice- versa). They satisfy [5]: 



N N N N 

Uk{x) K + h^Uxix) + hdu^K = 

N 

N N N N 

y_,j{x) k + kVj{x) - hd^jA = . 

N 



(2.36) 
(2.37) 



The expressions of the matrices J7,/, Vk in terms of P and Q can be found with the same notation in X but are 
essentially well known since |22| in a slightly different context and with different notation. 
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2.2 Deformations with respect to h 

The /i-deformations of the biorthogonal polynomials (derivations w.r.t h) are also represented by semi-infinite matrices 
of finite band size which we can write as follows: 

h^^n^ = H^^l, h^dn^^H^. (2.38) 

oo oo oo oo 

The relationship with the matrices P and Q is easy to derive but since we could not find it in the literature (except 
in special cases [201) we give it hereafter. From J ^'oo(a:)<&oo*(2/) = 1 we promptly obtain 



Moreover we have 



H ~QP + H^ ^0 . (2.39) 

{H-Vi{Q))^=0=(^H-V2iP'))^ (2.40) 

lower ^PPer 

H -Vi{Q) + {Vi{Q)~QP + V2{P))+H* -V2{P)^0 (2.41) 

{H - yi(Q))o = (^H- V^iP'))^ - iV,{Q) - QP + V2{P))o (2.42) 

iH-Vi{Q))_ = {QP-Vi{Q)-V2iP))_ (2.43) 

(^H-V2{P'))_ = {P'Q' -Vi{Q')-V2{P'))_ (2.44) 



so that we have 



H = Vi{Q)+ + ^{QP + Vi{Q)-V2{P))o + {QP-V2{P))-=Vi{Q)+o + {QP)-o-V2{P)-f, (2.45) 

H =. V2{P%o + {P'Q')-o-Vi{Q*)-o ■ (2.46) 

It is also interesting to consider another derivation of (|2.45ll . from scaling properties. Consider the change of 
variables x i-^ j and y i~+ ^ , it gives 

{Xfir+^hn{{uK}, {vj}, h)6nm = J dxdyA"7r„ {uk}, {vj}, h) /i™a„, {vj}, e-"'*''^1'''""" (2.47) 
where 

Vi{x) X^iVi (I) , V2{y) := A/zl/2 (^^^ , h = X^h . (2.48) 
This immediately implies the identities 

7r„(a;; {wkA^-^m}, {«j//-'A}, X^in) = A"^„ {u^}, {vj}, ft) (2.49) 



a„(y; W'A^-«M},{«JAi 'A},AAift) -//V„ J, {u^}, {t;,;}, (2.50) 

hn ({ukA^-^a^}, {vjAi^-^A}, A/ift) = (A/i)"+ift„ ({?/k}, ft) (2.51) 

V-n (x, {w^A^^-^A*}, {w/M^^'^A}, A^ift) = A'V^a^-'t^ V« {uk}, {■".]}, ft) (2.52) 

(/-n (y, {Mif Ai-^/i}, {wjm'^-'A}, A/xft) = Ai"^ A-'^0„ {uk}, {vj}, nj (2.53) 

In infinitesimal form, from Xdw, , and ud,,, we obtain 
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1 + n 
1 + n 



+ Y] UKduK - y^iJ - '^)vjdvj + fidn ^{x) = 

^ — ^ ^ — ^ / oo 

K J / 



K 



n — 1 



+ X! ~ X!^"^ ~ ^)'V]dvj + $(y) = 



\ K J J 

where n := diag(0, 1, 2, 3, 4, . . .). Using eq. (|2.3HI we immediately obtain that 



2 



(2.55) 
(2.56) 
(2.57) 



H = ^^"^ - PQ - ^(i^ - l)ui^t/K - ^ vjK/ = -^^i!!±ll + ^ ua'C/k + ^( J - l)vjVj (2.58) 

J K J 



which, using 

(1223 E21I), is identical to (1233). 
We also get the identity 

- PQ = Y'^^-'^'' + Y^^'<^'< '^"^ ■ (2.60) 

.7 /f 

This last relation is nothing but the expression of the scaling properties of the biorthogonal polynomials and will prove 
itself very useful later on. 



2.3 Ladder recursion relations (AE's) 

We give here only the main statements without proofs and refer the reader to ref. ISI for further details. Let us 

N 

introduce the sequence of companion-like matrices a(x) and a(x) of size ^2 + 1 



a(a;) 

N 



N 

1 







a(a;) 



'•• 

1 

1 



liN) 

z-^o(Af) 
7(JV-1) 

-ai(jV+l) 
7(JV-1) 



TV > do 



•• 

1 





7(Ar-l) 



iV > 1, 



N 



and also the analogous sequence of matrices hiy) and h{y) of size di + 1 



AT 



N 





1 






-fei(A') 



7(Ar) 
7(JV-1) 
7(iV-l) 






1 

-Pi(N) iy-f3o{N)) 
j{N) 7(JV) 

1 



N>di 



•• 

1 



7(JV-i) U U U 



, TV > 1. 



(2.61) 



(2.62) 



(2.63) 



(2.64) 
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The equations in H2.16|l imply the foUowing 

N ^ 

Lemma 2.1 The sequences of matrices a, b and a, b implement the shift N i-^ N +1 and N i-^ N —1 in the windows 

N N 

of quasi-polynomials and Fourier-Laplace transforms in the sense that 

N N+l jv 

a^^ix) = ^ ix) , Mx) = ^(x)& , (2.65) 

N N N+l 

N N+l N 

h^y) = $ (j/) , My) = * (j/)b . (2.66) 



N N N+l 



and in general 



N N+j N+'i-l N 

$ = a •••a^, a •••a, (2.67) 

N+j N+j-1 NN 

N N+j N+j -1 N 

$ = b •••b$, 3^=* b •b, (2.68) 

N+j N+j-1 NN 

N N — — — 

where ^ , $ here denotes a window in any of the Fourier-Laplace transforms defined in eqs. \2.1'^) . 

It was proven in that the differential-deformation-recursion relations H2.25l[TMIl2.65|l are compatible and admit 
a common fundamental solution. 

2.4 Folding 

Important remark. From this point on we focus on the wave- functions ipnix) and their dual '/'^^(a^), but everything 
being said can be immediately extended to the wave- functions (j)n(y),ijj^{y) by interchanging the roles of the spectral 
parameters, the potentials and the matrices P and Q. 

We recall that the notion of "folding" is the following: to express any quasipolynomial ipn{x) as a linear combination 
of (^2 -1- 1 fixed consecutive quasipolynomials V'JV-j , j — 0, ■ ■ ■ d2 with polynomial coefficients. We now provide a way 
of computing the folding which is different from (but equivalent to) the one used in W where the ladder matrices were 
used instead. 

Let N (i.e. a window) be fixed; we seek to describe the folding of the infinite wave- vector v[/(x) onto the window 

oo 

at TV by means of a single matrix J^{x) of size oo x {d2 + 1) and with polynomial entries such that: 

N 

N 

Vn, Mx) = V Tn,k{x)i^k{x) . (2.69) 

^ — ' N 

k=N-d2 



In fact it is more convenient to think of ^ as a oo x oo matrix with only a vertical band of width d2 + I of nonzero 

N 

entries (with column index in the range N — d2, . . . , N). In order to describe as explicitly as possible the matrix we 

N 

first introduce a convenient notation to express the band-matrices F, Q: 

Q = 7A + 2 ajK'' , P = A*7 + ^ K'Pj (2.70) 
i=o j=o 

where A is the shift matrix defined in (|1.7|) , and aj , /3j , 7 here above denote the diagonal matrices 

aj := diag(aj(0), . . . ,aj(n), . . .) , (3^ := diag(/3j(0), . . . , /3j(n), . . .) , 7 := diag(7(0), . . . , 7(71), . . .) . (2.71) 

It is known that the existence of biorthogonal polynomials satisfying H1.5|l is equivalent to saying that 7(n) 7^ 0, Vn. 
Therefore we assume that 7 is invertible. 
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Eqs. ifT^ imply that 

ad,A*'^ =t;,,+i(A*7)''^ (2.72) 

or, more transparently, 

{n) = Vd2+i-f{n - 1) • • • 7(71 - ^2) . (2.73) 
Let us now define the following semi-infinite matrices: 

A :^l ^(7-iA)^'^(g-x) , B:=l-n-(AV')(Q-a;) . (2.74) 

Vd2 + 1 

We note that A is strictly upper triangular while B is strictly lower triangular. 

When acting on the semi-infinite wave-vector ^'(a;) by definition we have (Q — x) \p = so that: 



* = A^', ^^[B + U]^ (2.75) 

oc 00 00 y J 00 

Notice that 1 — A and 1 ~ B are invertible because they are upper (resp. lower) triangular with ones on the main 
diagonal. Notice also that the matrix Q — x has a left inverse'*: 

^— (l-A)-i(7-iA)'^^(g-x) = l (2.76) 

and a right inverse (write 1 — i? = n + L*7^*((5 — x), and multiply by 7L on the left and by (1 — B)~^L*j~^ on the 



right): 

(Q-x)(l-S)-i(AS-') = l (2.77) 
but they do not coincide (one is upper and the other is lower triangular). 

2.4.1 Upper folding 

We start our construction of the folding matrix by first looking at the "upper" folding, i.e. the part of T with first 

N N 

index greater than N. 

We introduce the following notations 

N M M N M 

n:=i-n, n:=n-n = n-n = nn, (2.78) 

N N N M N 

and we remark the following formulas which hold since B is strictly lower triangular 

N N N N 

nsn = 0; uBn = Bn. (2.79) 

N 

With these formulas in mind we compute 

N N N / N\ N N N N N /N N N ' 

n* = nB*==nsn + nU' = nsn* + nsn* = nsn^' + nsnsn* + nBn* 

C30 00 \n J OG N 00 00 N 00 \ N 00 00 

N N N N N „ N 

= (l + B)nSn* + n-BnBn* = (l + B)nSn^' +B^ n * = [iterating r-times] = 

N OQ 00 N 00 OQ 

= Vs^nBn^+s'^+'n* . (2.8O) 

— ^ N 00 00 

j=0 



■^The fact that Q — x has a left inverse (Q — x)j^ ^ is not in contradiction with (Q — x) VI/oo = 0. 

Indeed, (^{Q — x)£^{Q — x)^ {Q ~ {{Q ~ ^) ^oo) because both LHS and RHS are infinite sums, and the order of summation 

cannot bo exchanged. 
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Given that B is lower triangular (in a sloppy sense it is nilpotent) the remainder in the formula above is a vector 
whose nonzero components are only those with index greater than N + r + l. Since in our problems we are only folding 
finite band matrices and we will never need arbitrarily high indices, we can disregard the remainder and send r —f oo 
(i.e. we take an inductive limit). The result is then 



n* = (1 - B)-^nBn^ 

oo N oo 



(2.81) 



Finally note that, since the matrix B has only d2 bands below the main diagonal, we have 

. 



N N 

nBn n 

N N-d2 



In other words the expression in ea. H2.81|l contains only the quasipolynomials ?/'JV-d2J- 
achieved the first part of our computation 



(2.82) 

,ipN- Therefore we have 



N , N 

n = (i-s)-i nB n 



N 



N 



(2.83) 



which can be simplified further as follows 



JV 

N 



(i-B)-i n B n - 

JV 

-(l-B)-i n(A*7"i)(Q-a:) n = 

N 

-(1 - B)-^ (AV') "n' (Q - x) n = 

JV 

JV-1 

-{i-B)-'{A'j-'){Q-x) n n-(i-s)-^(AV') 

JV 



Q~x, n 

JV-1 



n 

N 



JV-1 



= -(1 - B)-' 1 1 - s - n ) -(1 - B)-' (AS"'^ 



JV 



Q, n 

JV-1 



^-1 

- n -(i-B)-i(AS"') 



Q, n 

JV-1 



(2.84) 

Concluding this part we have 



JV JV-1 , , , 

UT = - U -(l-B)-i(AV ) 

JV JV 



Q, n 

JV-1 



(2.85) 



2.4.2 Lov^rer folding 

Similarly to the case before we have now 



JV-da-l 

n A u 

JV-d2-l 



; U A U 

JV-d2-l JV-d2-l 



A n 

JV-d2-l 



(2.86) 



and hence we promptly obtain along the same lines (we set for convenience M ~ N ~ d2 ~ 1) 



M 

JVf oo J\/f oo Af oo M M oo 

M ( M 

= nAn^' + nAnAn* + nAn* 

M oo M \M oo M M oo 

M „ 

= (l + A)nAn* n * = [iterating r-times] 



M 



M oo 



M 



V n A n * +A''+^ n 

j=0 



AI CO 



(2.87) 
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For the same reason used above we disregard the remainder term and thus we obtain 



M oo M oo 

Once more this expression contains only the quasipolynomials within the window at N. Therefore: 

iiT = {i- A) ^ n An 

M N M 

Once and again this expression can be simphficd further as fohows 

njF = {I- A) ' u An = 

M N M 

= — —{i-A)-^ii{n'^Kf-{Q-x)n^ 

Vch + l M 

= — —{i-A)-^{j-^Af- n{Q-x)n = 

Vd2 + 1 N-1 

= — —ii-A)-'i-f~'Af^Q-x) n n 

Vd2 + 1 N-1 

+ ^ (1 - A)-' ii-'Af^ [Q-x, n ] = 

Vd2 + 1 N-1 

= -{i-A)-^{i~A) n 

+^{l-A)-^{^-^Af^[Q, n ] = 

Vd2 + 1 N-l 

= - n +^{i-A)-^{^-^Ar^[Q, n]. 

N-l Vd2 + 1 N-1 



(2.88) 



(2.89) 



(2.90) 



FinaUy we remark that ij for N ~ d2 < i, j < N is the c?2 + 1 identity matrix on the selected window so that, putting 

_ N 

eqs. H2.85I I2.90|l together we have: 



N 



N M 
M N N N 

-(l-i?)-i(AV') 



Q, n 

N-1 



1 



(1- A)-i(7-iA)'^^ 



Vd2 + 1 



Q, n 



(2.91) 

We can summarize the above computations into the single formula 




(2.92) 

Using this folding operator we can immediately find the folded version of any other operator, viz, if L is any semi-infinite 
matrix with finite band size, its folded counterpart is: 



M 

L — > C[x) ^ ULTn 

N 



(2.93) 



Note that -strictly speaking- formula (|2.93|) defines a semiinfinite matrix with only a nonzero square block of size 
(i2 + 1 located on the main diagonal 

C{x)^j=Q ,i,j<^N~d2,...,N . (2.94) 

By abuse of notation we will denote by the same symbol the (o?2 + 1) x (^2 + 1) matrix corresponding to the nonzero 
block. In this fashion we will not distinguish between tracing over the finite or the semiinfinite matrix (with some 
advantages which will appear later). 
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2.5 Dual Folding 

For completeness we report the formulas for the folding of the dual wave vector $ onto the window iV — 1 < n < 

oo 

N + d2 — 1. The computation is completely parallel to the above using now the matrices (remember that the dual 
wave vector is a row- vector and hence we must act on the right) 



B = l- {Q~x)L*-f- 



-lj\d2 



A^i- n -{Q-x 



Vd2 + 1 



(2.95) 



With this in mind it is straightforward to obtain 



N-2 , N-2 

$ n - $ n B n ii-B)-^ = - n 

ooN-2 oo N-2 N-1 

$ n n A n (i-A) 

oo ocN+d2 — l 



N N-2 

£ = $ n B n (i-B)-i = 

oo N~2 



Q, n 



?, n 

Ar-1_ 

n 



L*7~^(l-B)-i 
1 



-lV'(i-^)-i + 



Q, n 

1 



'Vd2 + l 



(2.96) 
(2.97) 



3 Folded deformation matrices Uk{x] 



N 



From eq. H2.85|l and using the folding relation 12.9311 we obtain the desired formula for the deformation matrices Uk 



(the folded version of -j^Q+q), indeed (recall M = N — 6,2 — 1) 



N 



KUk 

N 



M „ M „ M 

n Q%Tn = n qJo n t= 

M „ Af M „ N 

n Q+o n n n J^jv = 

N N N 

nQjon-nQ"^ n -n g^(i-B)"'(AV ) 



N N N 



N N 



Q, n 



M M M „ N-1 

n Q+o n - n g_o n 

-n(Q^-a;^")(l-B)-i(A*7"i) 



M 



-x"" n(i-B)-i(A*7-i) 

Af 



Q, n 

JV-l 



Q, n 

AT-l 



uQ% n - n g^o n - n - 

AT Af-^l Af AT Af ' 



K 



Q, n 

AT-l 



AT-l 

^ n 



A'/ „ M M „ AT-l Af 

n 0^0 n -iiQ% n -nM^i^(x) 

AT Af-1 Af N N 



Q, n 

N-l 



„K 



N~l 

n 

AT 



M „ Af Af 

nQ% n -nw^if(x) 

AT AT-l AT 



Q, n 

Af-l 



^ 1 \ 

2^ ~ t^Qn.n I n , 

z In 



(3.98) 



where we have defined 



Wk{x) := 



(3.99) 



This formula generalizes similar formulas for the one-matrix model appeared in ^ which allowed us to perform explicit 
computations linking the isomonodromic tau function and the partition function of the model (see Conclusion) . Trace 
Formulas 

We compute the trace of the deformation matrices by noticing that the last term in H3.98(l has zero trace because 
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Wk{x) commutes with Q and because of the cychcity of the trace 



M 

Tr UWk{x) 

N 



Q, n 



Q, n 



M 

n = TrWKix) 

N 





Q, n 



= 



where Tr means the trace of the (finite rank) semi-infinite matrix. Therefore 



1 



K tr UKix) ^x"^ - ^Q%j^ . 

N I ^ — ' I 

m— 1 



0.1 -. 



K 

N-m,N-m 



(3.100) 



(3.101) 



Now, the deformation equations imply the deformation equation for the normafization coefficients of the biorthogonal 
polynomials, 

(3.102) 



which implies promptly 
trUKix) 



N 



In 



K 



-X 



(3.103) 



where we have used the fact that [TB| 

N-l 

Zn^CnYI hj , (3.104) 

3=0 

with Cn dependent only on N but independent on the V^'s. 

4 Differential system Di{x) 

N 

From the formulas for the deformation equations we can derive the formula for the differential equation. Recall that 
P is of finite band size and that we have: 



K=0 ^ 

Using lfTn|l and ^^H^ we get: 

di di 

F = L*7 + uil - J2 Kuk+iUk + 2 H ^k+iQo ■ 



(4.105) 



(4.106) 



K=l 



Since D (x) is nothing but the folded version of P, and the folding is linear, we have: 

Ni 



^ di 



Di{x) ^ IlL 'fJ^N + uild2+i - y uk+iKUk{x) + - y uk+i1 

N N f-^, N 2 Nn 



where 7 and <? are the diagonal matrices: 

N Nq 



7 =diag(7(^-d2-l),...,7(n-l)) , q'^ = diag (Q 



N 



Nr 



N-d2,N-d2^ ■ ■ • ' ^N,N) 



We leave to the reader to check that the first term is 

/liM) 



M . M . M M . 

nL*7jr = nL*7n + nL*7nj^^ 

N N N N N MN 



/ ad3-i(jV-l) 
1 



V 



7(7V-1), 



a:-Q:o(Af-l) 7(jV-l 
ad2{N-l) aa2(N-l) 





(4.107) 
(4.108) 



/ 

(4.109) 
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where the second matrix in the product is nothing but the inverse of the ladder matrix a (x) . Using the expression 

N-l 

for the folded matrices Uk and with simple manipulations one obtains {M = N — d2 — I). 



N 

N M MM 

Di(x) = n^i'(Q)> n - liWix) 

N - N-1 N 



Q, n 

JV-l 



N-l M . 

+ vi{x) n +nL*7J^Ar 

N N 



where W{x) is the semi- infinite matrix (polynomial in x and Q): 



(4.110) 



(4.111) 



In a totally explicit fashion we have the expression of the matrix Di(x) 



Di{x) 

N 



N 



/V{{Q)N-d2,N-d2 ■■■ V({Q)N-d2,N-l \ 



+ 



/7(M) 



V 



VI{Q)n-i,n-i 

V{{x)/ 

1 



7(iV-l), 



V 



/W{x)m+1,N-1 ■■■ W{x)M+l,N+d2-l' 



N 



V Wix)N,N- 



W{x)N,N+d2~ 



which is the exact analogue of the very explicit formula appeared in jS]. 
Notice that (using and if fc > we have: 

VliQ)N-k,N-l = -PN-k,N-l = (ik-i{N - I) 



/ 



(4.112) 



(4.113) 



which allows to rewrite: 



Di{x) 

N 



/'Po{N-d2) ... /3d,_i(A^-l) \ 
j{N-d2) ': ': 



V 






7(M) 



ad.(M-l) 

/ W{x)N-d2,N-l 



MN-1) 

7(iV-l) V{{x)J 
/ad2MN-l) ... ao(^-l)-a; 7(A^ - 1) \ 
... 

... / 



V 



V W{x)n,n- 



W{x) 

N-d2,N+d2-l 
W{x)N,N+d2-l 



N 



Trace formula 



(4.114) 



Using the previous formulas for the deformation matrices we can promptly obtain the trace of the matrix Di. 

N 



d2 



trDiix) = Vlix) + V/3o(A^ ~ k) - ^{M) 

N 



k=l 



ad2^i{N~\) 
a<j,(iV-l) 



(4.115) 
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and using eq. H2.21|l we have: 



d2 



ad. 



k=l 



and 



therefore: 



d2-l 



d2 



fe=i V j=i 



tvDi{x)^Vi{x)- 

N Wda + l 



(4.116) 



(4.117) 



(4.118) 



Note the remarkable fact that the trace of Di does not depend on N . We also point out that this result cannot be 

N 

obtained directly from the formula contained in [5]. 



4.1 Folded deformation matrices Vj{x) 



N 



We first consider the deformations w.r.t. vj, < J < c?2 because they have simpler expression and allow to compute 
the deformation w.r.t. Wda+i by folding of H2.60|l . The semiinfinite deformation matrix is the finite-band size matrix 



Its folded counterpart as per eq. (|2.93|l is (for brevity in the formulas we set M := N — d2 — 1): 

M , M 

-JVj{x) = n(P^)-o^jv-n(P'')-o n^AT 



N 



N 



N 



N 

N N N M 

N N N M 



From the expression 12.74|l for A we have the relation 



M 

n An 

M 



^ n(7"'A*)''^(g-a;)n 



Vd2 + l M 



^ -(7-iA*)'^^ n {Q-x)n 



Vd-2 + l 



N-l 



Vd2 + l N-1 

= -il-A) n +^(7"'A*)" 

N-l Vd2 + 1 

Plugging H4.121|) into H4.120|) we obtain 



Vd2 + 1 

Q, n 



Q, n 



JV.j{x) 



N 



M , M M , M 

nPioH-nP'^ n 

N N N N-l Vd2 + 1 N 



np-^ii- A)-\-f^^Af^- 



Q, n 

N-l 



(4.119) 



(4.120) 



(4.121) 



(4.122) 



The last deformation matrix is obtained from eq. H2.60|l by also recalling that jj^^+i QJ^^ — xU^^^^ Tjy to be 



di + l 



(d2 + l)vd2+iVd2+i ^ xDi{x) + J vjVj + KukUk - Miag(7V - d 

N N -, — N TT^ N 



(4.123) 



J=i 



K=l 



4.1.1 Trace 



We will compute first the traces of the deformation matrices for J < d2, whereas the case J = c?2 + 1 will be obtained 
by tracing the identity H4.123|) . Hereafter the notation tr denotes the finite dimensional trace for the folded matrices 
(of dimension (^2 + 1) x (^2 + 1)) while Tr denotes the trace of a semiinfinite matrix. Bearing this in mind we compute 



1 / jN-d2-l\ ( jN-d2-l\ 1 / , 1 1 J 

JtTVj{x) = -TT [PJ n - Tr P'^ n + Tt [P'' {l-A)-'{-,-^Af- 

N 2 \ N J \ N-l } Vd2 + \ V 



Q, n 

N-l 
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Af-1 



■j=M 

Let us focus our attention on the last term. First, notice that 

1 



Q, n 

N-l 



(4.124) 



1 



Tr P'^ (1 - [V2'(P) - X, Hat^i] 



(4.125) 

where A is the foUowing strictly upper triangular matrix 

Indeed, for J < d2, the difference is the trace of a strictly upper triangular matrix. Then we have: 



TrP^(l-yl)-i(7-iA)'^^ 



V^iP)~x, n 



N-l 



^ TrP" 



(l-A)-l(7-lA)'^^y2'(P)-.T 



n 

N-l 



(4.126) 



(4.127) 



Let us denote by yj ~ yj{x), j ~ 1, . . . ,d2 the solutions of algebraic equation V^iy) ~ x and define the following 
polynomials in y: 

^2' (2/) - X 



y-vj 

We have in particular Sj{yj) = V2'{yj). 

We also define the following strictly upper triangular matrices: 

R, :=l-(7-^A) {P-y,) 

We are now going to use the following identity, whose proof can be found in appendix IXI 

d2 



(1 - A)-\^-'Ar^ = (1 - Rjr'i-'A 



Using H4.130|l we have: 



Vd2 + 1 
d2 



TrP^ 



{l-A)-^{j-'Ar^,V^iP)~x 



n 



^ Wi^) " "^'^'^ ^'^''^ ^ 



N-l 



- E Tr P' iVm -x)il- P,)- 7-A = 

=1 

E Tr P^ (1 - P,)-i 7-^A(P - y,)^fl^ U 



d2 



d2 



d2 



y 7777^ Tr P^ S*, (P) n ~Y —7^ Tr P'^ Sj (P) A* A n 



N-l 



(4.128) 



(4.129) 



(4.130) 
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E vkr^ ^^■(^) n(P - y,)(l - 7-'A n 



, = 1 ^2"(%) 



d2 



(4.131) 



N-l 



Notice that if > 2d2 the Hn-i in the last trace is irrelevant and therefore that trace is independent of N. Assuming 
thus N > 2c?2 we obtain can carry on with our computation 



— Tr P' (1 - A)-' il-'Af' [Q, Un-i] 



Vd2+1 

d2 



d2 1 d2 



'^^'^'^^^''^ n(P-.,)(i-P,)-7-A5,(p) 

- E TTTtJ-T Tr n(P - %)(1 - R,)-' l-'A (P'' - y/) 5,(P) 



d2 



d2 

E 



Trn(P-2;,)(l-P,)-S-^A5,(P) 



P'^ - yj 

' SAP) 



•^2 ,,J 



E Tr p^ 5,(P) n - Tr n(P - y,){i - R,)-' j-'A (5,(P) - s^y,)) 



d2 



d2 



Y^vi Tr (P - y,)(l - R,)-^ 7-^A n - J] Tr n(P - y,)^^^^^ 5,(P) 



Recalling that A n = we continue 





d2 



U^i'iVj) 



™ - E^TrP.^.,(P)n-E^Trn(P'^-,/)5,(P) 



<i2 J 



(i2 d2 

Ey/Trn^j:^/ 



With eq. H4.133fl we can finally state 




(4.132) 



(4.133) 



(4.134) 
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Using eq. (|4.134|l and tracing (|4.123|) one can obtain the following formulas for the traces 



tr V,/(a;) = In 
N 2 



N-j 



In 



dy— 777^^ > J=l,...,d2 



J Viiy) - X 



tr Vd^+i{x) = --<9^, In 



N 



In 



dy 



d2 + 1 V^{y) - X 



-h N-^ 



1 



2 / Vd2+i 



(4.135) 



4.2 Determinant of the fundamental system 

It was proven in that the following overdetermined set of differential equations for a (^2 + 1) x {d2 + 1) matrix v[f 

N 

is Frobenius compatible: 



-Hd^^ix) ^Di{x)^{x) 

N N N 


(4.136) 


hduK *(a;) = Uk{x) ^{x) 

N N N 


(4.137) 


-hd^, ib{x) = Vj{x)<^{x) 

' N N N 


(4.138) 


a{x) ^{x) = \I> (x) 

N N N+1 


(4.139) 



The biorthogonal polynomials constitute one column of \I>(x), while the other d2 can be also explicitly obtained as in 

N 



From (|4.135|l it follows that the dependence on the coefficients of V2 of the determinant of the fundamental system 



IS 



d2 + l / ^ d2 + l 

Mvlndet ( *(a;)j := - X] Indet [f^{x)) = ^ dwjtr ( Vj(a 



- hdv In 



N-j 



IN 



d2 + l ^ 

^dvj j> dy 

J=i 



.1=1 



N 



y' n'jy) 

J V^{y) - X 



Now one can check directly that 



J V^iy) - X - '''' 



dy iV2{y) - xy) — 



Vi'{y) 



(4.140) 
(4.141) 

(4.142) 



Vi{y)-x, 

Therefore, as a consequence of eqs. ()4. 1181 IT. 1351 13. 1(1!^ we find that any matrix fundamental solution to the compatible 
system of equations satisfies the equations 



hd^lvL det*(x) = -V[{x) 



N 



Vd2 
Vd2 + 1 



Muln ( det*(:r) ) ^ d„ l-V,{x) + J In ( ^"4^^'' 



N 



hd-v In [ det "^{x) ) = —hdv 



In Vd2+i ^ 



-TV, 



IN 



+ fdy {V2{y)-xy) — 



V^'iy) 



V^{y) - X 



N+l N 



det * (x) ^{xy = det(a(a;)) = Vd2+i{-) 



d2 + l / hN~d2 



N 



IN 



and hence we have finally the complete formula 

{Vd2 + lf-'^ (_l)A^('i^ + l) W nl=l hN-m 



det 



*(x) 

N 



T-N 



exp 



-T[Vi{x)+idy{V2{y)^xy) — 



Vi'{y) 



Vi{y) - X 



(4.143) 
(4.144) 

(4.145) 
(4.146) 

(4.147) 
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Note that the contour integral just returns the sum of the critical values of ¥2(1/) — xy. 

We know from 5 that a joint solution of the dual system of overdetermined equations for the dual window denoted 

N 

by ^{x) has the property that 

N N 

A*(a;) = C , (4.148) 

N 

where C is an invertible matrix which does not depend on any of the potentials ot n or h or x and can be conveniently 
be normalized to unity (see 0| for an explicit construction). Therefore 



N AT 

det(*(x)) = det(A)" det(*(x))"^ 

N 



(4.149) 



5 Dual folded system 



N 



For completeness we add the formulas for the relevant folded operators for the dual window $ = . . . , 4>N+d2-i)- 

Since the steps are essentially the same we give only the results (here we set L = N + d2 — 1) 



6 



N 

KUk 


N- 

= n 

L 


1 

Q 


N 

JVj = 


= n 

L 


Pi 



N-l 



N~l 



n Q% n -[Q,nN-i]WK{x) n 



N-l 



N 



\ 



l{L) 



adiiL + 1) 

Conclusion 






Vo 



5« 



K 

N-1,N-1 



. N-1 







/ V{{x) 
jiN-1) (3oiN) 



iiN-l)\ 

aa{N) - X \ N 



ad2-i{L) / 



/ W{x)N~d2,N-l 

V W{x)n,n~i 



N-l 

n 

N 



W{x) 

N-d2,N+d2-l 
W{x)N,N+d2-l 



(5.150) 
(5.151) 

(5.152) 
(5.153) 



We have given the most explicit construction for the matrices describing the differential-deformation-difference folded 
system. This has allowed us to determine explicitly the determinant of the fundamental solution of the system in 
terms of the partition function of the model and the two potentials. The final expression is quite simple in comparison 
with the complexity of the computation, expecially for the traces of the deformation matrices Vj. 

N 

It is our hope (in fact it is our plan) that these formula be used to relate explicitly the partition functio of the two- 
matrix model to a (suitably defined) isomonodromic tau-function. Indeed the system (|4.13(jl4. 13714. 13814. 1391) can 
(and should) be regarded as a monodromy-preserving set of infinitesimal (|4. 1371 B. 138(1 and finite H4.139|l deformation 
equation for the ODE (|4.13t)|) which has an irregular, degenerate singularity at infinity. As such it is envisionable that 
one can define the associated notion of tau-function although there are technical difficulties due to the degeneracy at 
infinity. In the similar case of the ODE-|-PDE-|-AE for the orthogonal polynomials in the one-matrix model a similar 
approach has given very satisfactory results |51 |U] . 
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A Proof of formula M.l^Oh 

We need to prove formula H4. 130(1 



(1 - A)-\r^KY^ = ^TTTZTT (1 - Rir'r'^ , (A.154) 
which can be written more transparently as 

We remark for the better understanding of the reader that if the matrices A, (V^iP) — x) were invertible this would 
amount simply to the partial fraction expansion of 1/(^2 (-P) — x). Multiplying both sides by the invertible matrix 
(1 — A) on the left and recalling the definition 14.126|l of A 



(v4.156) 



-y,)(l-i?,)-V'A 
+A*77-'A)(P-2;,)(l-i?,)-i7-iA 



V'(P) — X 

Noticing now that {'j^^AY^— n = we continue 



d2 



1 A^^^2 ( -1 _ TT I — ('^,-1 A^'i2 



(7"^A)''^ ( 1 - n ) = (A.157) 



where we have used that (7 A) ^ n = 0. The identity is thus proved. 
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